Asymptotic Scaling of the Diffusion Coefficient of Fluctuating "Pulled" Fronts 
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We present a (heuristic) theoretical derivation for the scaling of the diffusion coefficient Df for 
fluctuating "pulled" fronts. In agreement with earlier numerical simulations, we find that as A'^ oo, 
Df approaches zero as 1/ln^ A'^, where A'^ is the average number of particles per correlation volume 
in the stable phase of the front. This behaviour of Df stems from the shape fluctuations at the very 
tip of the front, and is independent of the microscopic model. 
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Interest in the effect of fluctuations on propagating 
fronts has revived in recent years with the understanding 
that fronts which in the deterministic mean-field limit 
are of the so-called pulled type are surprisingly sensitive 
to fluctuation and discrete-particle effects. Pulled fronts 
are those that propagate into a linearly unstable state, 
and whose asymptotic front speed t^as is simply the linear 
spreading speed v* of infinitesimal perturbations around 
the linearly unstable state 0, S Si- The propagation 
mechanism of such fronts is that they are being "pulled 
along" by the growth and spreading of small perturba- 
tions around the linearly unstable state. It was first es- 
tablished by Brunet and Derrida ^] and later confirmed 
in a variety of stochastic front equations 0, IE IE S 
that when pulled fronts are realized by stochastic moves 
of discrete particles on a lattice, such that the average 
number of particles per lattice site or correlation volume 
in the saturation phase of the front is N , v^a approaches 
V* from below extremely slowly: the convergence to v* 
scales only as l/lr? N when N ^ oo, with a known 
prefactor that depends on the model under considera- 
tion. The reason why a true pulled front is so sensitive 
to finite particle cutoff effects is the very fact that there is 
essentially no growth below the cutoff level of one "quan- 
tum" of particle. As a result, these discrete particle front 
realizations are actually weakly pushed |EiEjO|- To re- 
mind ourselves that they converge to pulled fronts in the 
limit N —* oo, here we will refer to them as fluctuating 
"pulled" fronts. 

The second important feature of fluctuating "pulled" 
fronts, namely the diffusive wandering of the front it- 
self around its average position as a result of stochas- 
ticity in the microscopic dynamics, still remains very 
poorly understood. Of particular interest is the ques- 
tion how the front diffusion coefficient Df vanishes as 
N —^ 00. This scaling is particularly difficult to study nu- 
merically. Nevertheless, using a clever algorithm to take 
N as large as 10^^", Brunet and Derrida ^ presented 
convincing numerical evidence that for the model they 
studied, Df scales as ~ I/ln"^ TV. Moreover, a "simplified 
model" , where the fluctuations are randomly generated 
only on the instantaneous foremost occupied lattice site 
(i.f.o.l.s.), was found to exhibit the same asymp- 
totic scaling of Df as the full stochastic model . 



The microscopic dynamics of Brunet and Derrida's 
(discrete-time) model ,Sj closely resembles that of the 
(continuous-time) "clock model" 0. In both models, one 
considers a set of N particles with integral "readings" k = 
0, 1,2,.. .. The number of particles with a certain read- 
ing k at time t is nk{t). With (j)k{t) = J2'k'=k''^k'{t)/N, 
in a computer simulation, the speed and the diffusion 
coefficient of the front in both models are measured by 
tracking the position of the centre of mass of the particle 
distribution (henceforth referred to as the centre of mass 
of the front itself) in individual realizations. For the par- 
ticle distribution {nk{t)}r of realization r at time t, in 
both models, the centre of mass of the front is located at 



(1) 



From there onwards, one defines Df ^ limj-^oo d{[Sr(t-{- 
T) - S,{t) - VNTf)/dT and vn = {Sr{t)). The angular 
brackets denote first an average over all possible updating 
random number sequences for a given front realization 
after time t 3> 1, and then a further average over the 
ensemble of (initial) front realizations at time t. 

To the best of our knowledge, at present, there ex- 
ists no analytical derivation of the 1 / In^ N asymptotic 
scaling of Df for fluctuating "pulled" fronts. Our aim 
here is to provide a derivation for the stochastic Fishcr- 
Kolmogorov-Petrovsky-Piscunov (sFKPP) equation 



dt(f> ^ dl(f> + (f> -c^^ + ^^^^7j{x,t)/VN, (2) 

where the stochastic term cx 77(2;, t) is interpreted in the 
Ito sense with the two following conditions: 

{•nix, t))^ = 0, (77(x, t)Tj{x' ,t')),, = 5[x - x') S{t - t'). (3) 

We then argue that the same scaling of Df holds for 
the clock model |3 and for the microscopic model that 
Brunet and Derrida considered for their simulation — 
in other words, the asymptotic scaling of Df is indepen- 
dent of the microscopic model and is a generic property 
of fluctuating "pulled" fronts. The full flavour of this 
subtlety-riddled derivation appears elsewhere (Sees. 2.5 
and 4.2, Ref. ,12j). Here we focus only on the main 
points and the main results. 
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From this perspective, it is therefore important to 
first summarize the Langevin-type field-theoretical ap- 
proaches for general reaction-diffusion systems [T^ ll^ 
Starting with the Ito stochastic differential equa- 
tion 



M = dl4> + /(0) + e^l'^R{x, t) , 



(4) 



with e" <C 1 and R(x,t) = g(<j))ri{x,t), in these ap- 
proaches, one writes the corresponding front solution as 

t) = [x - V^st - X{t)] + 5(l)[x - Ve.st - X{t)] to 

separate the systematic and the fluctuating part of the 
front. The systematic part satisfies the determin- 

istic equation -Was^^ (C) = d^(l>'^°HO + As 
for the fiuctuating part, the idea behind writing (j){x, t) in 
the above form is to separate the fluctuations in the front 
at two different time scales. Of these, the long time scale 
fluctuations are coded in the random wandering X{t) of 
the Goldstone mode $G,fl(C) = d(l)'-°\^) / of the front 
around its uniformly translating position x — Vg^t. On the 
other hand, the short time scale fluctuations manifest 
themselves through the fluctuations in the front shape 
around its instantaneous position x — v^st ~ X{t). 

In this form, the Goldstone mode is in fact the 
right eigenvector of the linear stability operator = 
9| -I- Vasd^ + ^0/(0)10=0(0) with zero eigenvalue and 
its instantaneous position is defined by requiring that 
it be orthogonal to the front shape fluctuations, i.e., 
I^d4<^G,L{0m - Xit),t] = at all times lH M 

Qllll- Here, <I>g,l(0 = e'''"'^$G,fl(C) is the left eigenvec- 
tor of Cv^^ corresponding to eigenvalue zero 0|. When 
this orthogonality condition is used in linearized Eq. 
(@J, one finds that the instantaneous speed of the Gold- 
stone mode fluctuates (at long time scales) around Was by 

Cllillillil 



X{t) 



(5) 



Simultaneously, at short time scales, the shape fluc- 
tuations of the front around cj)^^'^ at the instantaneous 
position of the Goldstone mode are analyzed by defin- 
ing the mutually orthonormal shape fluctuation modes 
{^'m(C)} in the eigenspace of non-zero eigenvalues of the 
linear stability operator as 0, IT3 | 



H{^,t) = ^c„(t)vl/„,fl(0. 



(6) 



The Cm{tys are then easily seen to satisfy [l^lT^ 

/"OO 

Cm(0 = - C^{t) + £-1/2 / *„,i(C)i?(e, t), (7) 



where 4'm,_R(^) and 5'm,L(C) ^-re respectively the right 
and the left eigenvectors of £^^^ with eigenvalue ^ 0. 

If the front position is defined by the position of its 
Goldstone mode, then the phenomenon of front diffusion 
arises due to the random speed fluctuations X{t) of the 
Goldstone mode around Uas- The diffusion coefficient of 
the Goldstone mode (which is identified with the diffusion 
coefficient of the front itself in this case) is then easily 
defined via the Green-Kubo relation |3 0, Q, as 



-^G = 2 — -r T2 — ■ y^i 



The angular brackets with a subscript t denote an average 
over an ensemble of (initial) front realizations at time t. 

On the other hand, if the front position is defined by 
the position of its centre of mass (Q, then in continuum 
space, we have S{t) = ^ _^dx (j){x,t), and for reaction- 
diffusion fronts Q that satisfy <t>''^\(,) ^ 1 for ^ ^ — oo, 
0(°HC) ^ for C ^ cx), and (50(0 ^ for ^ ^ ±oo [12 



S{t) 



df{4>) 



0(0) (4) 



diR{^,t). (9) 



f,i(t) 



b2{t) 



Thereafter [having used Eqs. H3I4|) . © and the solution 
of 0], in the expansion of the product S(t)S{t -I- t') of 
Eq. ()10|l . bi{t)b2{t + t') is seen not to contribute to 
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while d'^P ^ d'^P and D'^f , the respective contributions 
of b2{t)b2{t + t'), b2[t)bi{t + t') and bi{t)bi{t + t') to 



t+T 

dt' {{S{t)S{t + t'))r,)u 



l(3) 



(10) 
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Trn{Crn{t) Cm'{t))t 



diig^mMt 
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771, R 



(C) 



0(0) (5) 



d^ 



, Sf{<f>) 



*m',i?(r) 



0(O)(^'^ 



(11) 

(12) 
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In general, in this formalism, g{(p) cannot be replaced by ^[(/j^*^-^] — after all, the value of (j) at any given time 
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depends on its precise evolution history (i.e., noise realization chosen to evolve the front) at earlier times. However, 
in case the nonlinearities in /(0) make (j)^'^'^ pushed, i.e., for fluctuating pushed fronts with internal noise, the de- 
cay times of the shape fluctuation modes are of the same order as the time scale set by 1/was- In that case, the 
dependence of g{(j)) at any time on the precise evolution history of the front at earlier times can be neglected to 
replace g{(j)) by g[(f)^^^]. This replacement makes the stochastic (noise) term in Eq. ^ additive in nature. More- 
over, it converts Eq. to that of a Brownian particle in a fluid with viscosity r,^^, where the fluctuating force is 
given by e^^^ (i^ ^m,L('C)ff['/'^"''(C)]^/i implying that one can then use fluctuation-dissipation theorem to evaluate 



to further obtain 12] 



D 



(3) 



I.e., Df = — 



(13) 



We will see later that for fluctuating "pulled" front in 
the sFKPP equation too, Df and Dq are not the same. 
Before we delve deeper into the sFKPP equation, here we 
take a short digression to mention that similar situation 
occurs for gas mixtures (see for example. Chap. 11.2 of 
Ref. 0]). Therein, the expression (and the value) of 
the diffusion coefhcient depends on its precise definition, 
but these different expressions of the diffusion coefficient 
are (quite non-trivially) related by means of Onsager re- 
lations for the diffusion coefficients. As for fronts too, 
it should not be surprising that the precise values of D/ 
and Dg are not the same. Indeed what is important to 
note is that conceptually they are two entirely different 
quantities. Whether they could be related by any clever 
means or not is left here for future investigation. 

For the fluctuating "pulled" front in the sFKPP equa- 
tion, we make the ansatz that the front solution can be 
decomposed to a (/)*-°'' that is nothing but Brunet and Der- 
rida's cutoff solution, and its corresponding shape fluctu- 
ation modes {^r?x}; i-c, at the (linearized) leading edge 
of the front, c/y^^^O In A^sin[go(^ - 6)]e-^'V7r and 

^miO - sin[g™(C-ei)]e-^'Vln'^'^ 110. Here A* = l, 
qj = (j -I- l)7r/lniV Vj > 0, ~ is the location of the 
left end of the leading edge where the nonlinear 0^ term 
in the sFKPP equation is non-negligible, and the cutoff 
is implemented at — In iV |3j ■ This ansatz is con- 
sistent with the compact support property of the front 
solution in the sFKPP equation 18], and it also satisfies 
the requirements for Eq. Q [i.e., we can safely use Eqs. 
(|llll2f) with this ansatz]. 

Then, the 1/ln^iV asymptotic scaling of Dq for the 
fluctuating "pulled" front in the sFKPP equation is ob- 
tained very easily ,12] . One simply has to notice that due 
to the exponential weight factor in the integrand of the 
numerator of Eq. ©, the numerator is practically de- 
termined within a distance of 0{1/X*) of ^o, where the 
integrand scales ~iV and cancels the 1/A^ prefactor. The 
denominator, on the other hand, asymptotically simply 
scales as In^ TV. 

As for Df oi the fluctuating "pulled" front in the 
sFKPP equation, D^^^ and D^p are straightaway seen 
to be of 0{1/N) from Eq. lfTT|l . while the evaluation of 



Dj: is no easy matter. Unlike pushed fronts, one cannot 

simply replace g{(f>) by g[4i^-^^ and arrive at Eq. — 
for fluctuating "pulled" fronts, fluctuation modes decay 
very slowly wi, 131/ and as a result, the front configuration 
at any time depends strongly on the precise noise real- 
ization that has been used to evolve it at earlier times! 
This means that there is no other way forward than to 
evaluate D^p in Eq. (|12|l . and this is done in two steps. 

At the first step, we argue that in Eq. ((T^ 6vm = 
s„/ln^/^ N with Sm ^ 0{1) ^ and obtain 

= \ E rJ-^{c„Mc^,{t)),. (14) 

This scaling of 5vm is obtained with the idea that for 
Brunet and Derrida's cutoff solution 0^°^ , the front speed 
/.!1'^C/[<?!'^°H0] =i 2 is of 0(1). Naturally, when the 
front shape (/){£,, t) deviates from by an amount 

^m,_R(0; which is always a factor ln~'^^^ TV weaker than 
(0 itself [notice the prefactors of (/)(°) {£,) and "ifmMO 
above], the contribution of the m-th shape fluctuation 
mode to the fluctuation in the front speed Svm has to be 
of ©(In-^/^ as well [13 ■ Furthermore, since there are 
©(In TV) number of shape fluctuation modes [^0, the 
sum over m and m' in Eq. (|14|) runs from 1 to In TV. 

At the second (and perhaps the trickiest) step, we de- 
termine the dependence of {cra{t) Cm' it)) t on InTV, and 
evaluate the sums in Eq. H14|l . To this end, notice that 
for a given reafization, Cm{t) is expressed [via Eq. ^] 
as c„i(i) — ^'m,L <5</>(^, i) [13 • However, the pres- 

ence of e^ ^ in 5'm,L(C) implies that Cm{t) for a given 
realization is practically determined from the fluctuation 
characteristics at the tip of the front, and therefore, we 
retain the integral only over the leading edge of the front: 

^™(^) = TTTTT; /^"^^ sin[g™(e - 6)] H{L t). (15) 
in TV Jji 

By virtue of {6<j){^,t))t = 0, Eq. (HSJ then yields 
{c„i{t))t = as it should, but in the absence of any 
statistics of the shape fluctuations of the front at time t, 
one cannot obtain an expression of {cm{t) Cm' {t)) t from 
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it. Moreover, since we cannot replace g{(t>) by g[(j)^'^^ for 
fluctuating "pulled" fronts, one cannot use fluctuation- 
dissipation theorem in Eq. to evaluate {cm{i) Cm' {t)) t 
either — there is no generic fluctuation- dissipation theo- 
rem available for multiplicative noise flkl . Nevertheless, 
we can still proceed with two approximations. The first 
one stems from the fact that although it is clear from 
Eq. ifTCjl that Cm{t) and Cm'{t) {'m 7^ m') are corre- 
lated in general [after all, for a given realization, all the 
Cm(t)'s are determined through the same S(j){t)], these 
fluctuation modes will have a finite correlation "length" , 
i.e., {cm{t) Cm>{t))t will be negligibly small [compared 

to \/ {(^mi^))t{Cm' (0)t] when |m — m'j exceeds a certain 
threshold a <ti IniV. Based on this anticipation, our ap- 
proximation is to choose a = for the extreme (and un- 
realistic) case to simplify the expression for D^p to (see 
later for the discussion on non-zero values of a) 



In JV 



D 



(3) _ 



1 

It. 



(16) 



Then the second approximation is that due to the pres- 
ence of the ^ in the integrand of Eq. IjlSI) . only the 
value of 6(j) within a distance ~ 1/ A* = 1 of the tip deter- 
mines Cra{t)- This is seen in the following manner: typi- 
cally the magnitude of (5(/)(^, t) is of order y^(l)[^,t)/N ; at 
the tip, e^'^o ^ N cancels (50(^,t) - in Eq. (O, 
but further behind, the factor of y^4>{^, t)/N can 

no longer be compensated by e'*' ^. We therefore use 



I si n[g^(e-a)]| _ gm _ 7r(m + l) 



(17) 



In Eq. lfT7|) . sin[(7m(^ — Ci)] l^^s been Taylor expanded 
around its value zero at ^0 — In N . Thereafter, with 
T,„ = In^ iV/[7r2{(m + if - 1}] [li, we obtain 



one expects in reality, Eqs. lfTT|l . ((TH) . (fTK|l and lfT7|l show 
that the 1 / In^ N asymptotic scaling of I?/ continues to 
hold, (ii) We have extensively used the left eigenvector of 
the stability operator for reaction-diffusion systems. 
For clock model 0, or for the model that Brunei and 
Derrida considered for their simulation 0, construction 
of the left eigenvector for the corresponding Cy._^^, and 
repeating the same exercise (as in here) are nontrivial. 
Nevertheless, since all the arguments for the derivation 
of the 1 / In"^ N asymptotic scaling of Df in the sFKPP 
equation are concentrated on the leading edge (or more 
precisely, at the very tip of the front) where the (fluctu- 
ating "pulled") front properties are model- independent, 
one expects to observe the same scaling for D f for these 
two models too. In view of this, the 1/ln^iV asymp- 
totic scaling of Df seems to be a generic property of fluc- 
tuating "pulled" fronts, independent of the microscopic 
model, (iii) In the clock model, one can only create non- 
local fluctuations in the front shape 0. The "collisions" 
between clocks are also non-local in nature 0, [T^ . In re- 
ality, however, these complications matter neither for the 
front speed, nor for Df — Brunei and Derrida showed 
0, in a simplified version of their original microscopic 
model (which closely resembles the clock model), that 
the localized fluctuations in (/) at the very tip of the front 
is all that is needed for the 1/ In'^ asymptotic scaling of 
Df. (iv) Dg and [through the scaling of the €,„(<) 's] Df 
are both determined only from the tip of the front. This 
is in perfect agreement with Brunei and Derrida's "sim- 
plified model" 01 ■ (v) Finally, the scalings of D f and Dq 
have been obtained by means of an ansatz [second para- 
graph below Eq. H13() ]. The integrity of the method used 
here to obtain these scalings can be tested by numerically 
obtaining the scaling of Dq and the scaling properties of 
{cm{t)cm'{t))t for the FKPP equation. It must also be 
noted that these numerical simulations involve extremely 
high values of N , and are notoriously difficult to perform. 



Df^D 



InJV 



(3) 1 

/ o 2^ 



{m + l)'^ s"^ 



[(m + 1)2 - 1] In"* N In^ N 



(18) 



We end this paper with five final observations: (i) Real- 
istically, a 7^ 0, but so long as a <C InA^, which is what 
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